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Abstract. We consider heavy-quark energy loss and p±— broadening in a strongly-coupled N = 4 Super 
Yang Mills (SYM) plasma, and the problem of finite-extend matter is addressed. When expressed in terms of 
the appropriate saturation momentum, one finds identical parametric forms for the energy loss in pQCD 
and SYM theory, while p±— broadening is radiation dominated in SYM theory and multiple scattering 
dominated in pQCD. 



43 

r-| 1 Introduction 



PACS. 11.25.Tq Gauge/string duality - 12.38.Mh Quark-gluon plasma 
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It is often argued that hard probes such as heavy quarks 
are understood well enough to provide clean measure- 
ments of the properties of the quark-gluon plasma formed 
in Au+Au collisions at RHIC. In particular, they could 
help determine whether the plasma is weakly or strongly 
coupled. Results on bulk observables like the elliptic flow 
or the shear viscosity already prompted claims that the 
plasma is strongly coupled. 

With hard probes, it is indeed unclear if the pertur- 
bative QCD (pQCD) approach can describe the suppres- 
sion of high— p± particles, in particular for heavy-quark 
production. High— p± electrons from charm and bottom 
mesons decays seem to indicate a similar suppression for 
light, charm and bottom quarks. By contrast in pQCD, 
the heavier the quark the weaker the suppression. 

This motivates to think about strongly coupled plas- 
mas. The tools to address the strong coupling dynamics 
in QCD are quite limited, however for the N = 4 Super- 
Yang-Mills (SYM) theory, the AdS/CFT correspondence 
is a powerful approach. The findinds for the strongly- 
coupled SYM plasma may provide insight for gauge theo- 
ries in general, and some aspects may even be universal, 
like the lower bound of the shear viscosity pQ. 

In this work we study the energy loss of a very ener- 
getic heavy quark propagating through the SYM plasma. 
For comparison, the pQCD results are recalled in Section 
[2j The trailing-string picture of heavy-quark energy loss 
in a static infinite-extend plasma is introduced in Section 
[31 and a partonic interpretation in terms of the satura- 
tion momentum is given which allows to infer the plasma 
length dependence of the energy loss when considering 
finite-extend matter. 
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2 Heavy-quark energy loss ... 
weakly-coupled QCD plasma 

In this section we consider a heavy quark of energy E 
and mass M, propagating through a weakly-coupled QCD 
plasma and losing energy at the rate —dE/dt due to the 
interaction with the medium. We recall the pQCD results 
as well as the underlying multiple scattering picture. 



2.1 The heavy quark wave function 

Let us start with the heavy quark wave function in QCD. 
At lowest order in ot s , quantum fluctuations consists of a 
single gluon (see FigfTJ), whose energy we denote u> and 
transverse momentum k^ . The virtuality of that fluctua- 
tion is measured by the coherence time, or lifetime, of the 
gluon 

(1) 



t c = uj/kj_ . 



Short-lived fluctuations are highly virtual while longer- 
lived fluctuations are more easily put on shell when the 
heavy quark interacts. The probability of the fluctuation 



is 



P = P(M = 0)[ 1 + 



2 fc 2 



r 



with P(M = 0)~a s N c . (2) 
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We have introduced the Lorentz factor of the heavy quark 
7 = E/M . (3) 

Compared to massless quarks, the fluctuations with lu > 
jk± are suppressed in the wave function. This means that 
when gluons are put on-shell, they are not radiated in a 
forward cone around the heavy quark. This suppression 
of the available phase space for radiation, the dead-cone 
effect, implies less energy loss for heavier quarks [2]. 
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Fig. 1. Virtual gluon fluctuation in the heavy quark wave 
function, with energy ui and transverse momentum k±. The 
coherence time ujjk\ measures the virtuality of the fluctuation. 



2.2 Medium-induced gluon radiation 

In pQCD, medium induced gluon radiation is due to mul- 
tiple scatterings of the virtual gluons [3] . If, while undergo- 
ing multiple scattering, the virtual gluons pick up enough 
transverse momentum to be put on shell, they become 
emitted radiation. The accumulated transverse momen- 
tum squared pickep up by a gluon of coherence time t c 
is 



P± = A* — = 1 t 



(4) 



where fj? is the average transverse momentum squared 
pickep up in each scattering, and I is the mean free path. 
These medium properties are involved through the ratio 
q = fi /l, this is the only relevant information about the 
medium. In terms of the temperature T, one has 



q 



a,T 6 



(5) 



and at RHIC temperature, the value is q ~ 1 GeV 2 /fm 
(however 5 — 10 GeV 2 /fm seems favored by RHIC data). 

Since only the fluctuations which pick up enough trans- 
verse momentum are freed (k± < p±), the limiting value 
can be obtained by equating k\ with p\ 



qw/k\ 



kx < (qoj) 



1/4 = 



<2s(w) 



(6) 



The picture is that highly virtual fluctuations with k± > 
Q s do not have time to pick up enough p± to be freed, 
while the longer-lived ones with k± < Q s do. That trans- 
verse momentum Q s which controls which gluons are freed 
and which are not is called the saturation scale. 



2.3 Heavy-quark energy loss 

When applying this picture to heavy quarks, one sees that 
due to the dead cone effect, the maximum energy a radi- 
ated gluon can have is w = jk± = "fQ s , and its coherence 
time is t c — j/Q s . This, combined with the probability 
a s N c to have the gluon fluctuation in the wave function, 
allows to estimate the heavy-quark energy loss: 



dE 

~dt 



oc a<,N, 



= a s N c Q 2 s 



(7) 



The saturation momentum in this formula is that of the 
fluctuation which dominates the energy loss, and should 



be evaluated at u> = "fQ s - Using Q 2 S = \fquj one can ex- 
press that value of the saturation scale in terms of T and 
7 only: 

Qs = m l/3 . (8) 

In the following, that particular value is what we mean by 
Q s . We shall also use t c to denote the coherence time of 
the dominant gluonic fluctuation: 



7 2/3/£l/3 



2.4 The case of finite-extend matter 



(9) 



When the heavy quark is propagating through a plasma of 
finite length L, the discussion above has to be modified [4]. 
If L > t c , the gluons which dominated the energy loss in 
the infinite matter case have the time to form before the 
heavy quark exits the plasma, and therefore the matter 
is effectively of infinite extend and the above results are 
unchanged. If L < t c , only shorter-lived fluctuations, with 
less energy uj < Lk\ can contribute to the energy loss. 
The dominant ones have a coherence time uj/k\ — L, and 
the transverse momentum squared they picked up is qL. 
This defines the saturation scale in the finite matter case: 



Qt = qL 



(10) 



All gluons with k± < Q s are freed and the maximum 
energy of a radiated gluon is uj = Lk\ = LQ 2 S which gives 



— cx a s N c -^ = a s N c Ql 



(11) 



This formula is the same as before except that the satu- 
ration scale is different. In the finite matter case, Q s does 
not depend on E/M anymore, but on L instead. 

Note that the saturation momentum discussed here is 
not the traditional saturation momentum characterizing 
the small— x part of a hadronic wave function. qL has 
been denoted Q 2 S because when writting qL in terms of the 
gluon density per unit of transverse area in the plasma, 
one finds the same expression than when writting the sat- 
uration scale in terms of the gluon density per unit of 
transverse area in a hadron. 

2.5 Heavy-quark broadening 

From the above discussion, the radiative p±— broadening 
of the heavy quark is easily estimated: 



dpi 
dt 



ix. a s N c 



dQj 
dt 



a,N r 



(12) 



where t = t c in the infinite matter case and t = L in the 
finite matter case. The result is the same in both cases be- 
cause the broadening is due to multiple scatterings which 
give local p± kicks to the fluctuations. However radiative 
p±— broadening is irrelevant in pQCD (we mention it for 
future comparisons) with a s <C 1, because what is domi- 
nant is p±— broadening due to multiple scattering of the 
heavy quark itself: 



dt dt 



(13) 
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3 Heavy-quark energy loss in a 
strongly-coupled SYM plasma 

In this section we compute the rate of energy loss of the 
heavy quark —dE/dt in the strong coupling regime. To do 
so, we consider a N — 4 SYM plasma instead of a QCD 
plasma. The field content of this theory is 1 gauge field, 

4 fermions and 6 scalars, all in the adjoint representation 
of the gauge group. Using the AdS/CFT correspondence, 
the quantum dynamics of this theory at strong coupling 
can be obtained by classical gravity calculations [SB?]. 



3.1 The AdS/CFT correspondence 

We consider the large N c , small gauge coupling gyu limit 

N c -> oo , Qym -> , A = 9ym n c finite, (14) 

where the 't Hooft coupling A controls the theory. Then 
strong couping means A >> 1. In this regime the equiv- 
alent string theory in AdS$ space is weakly coupled and 
weakly curved: 



9ym and A > 1 <^ R > l s 



(15) 



where g s is the string coupling, l s is the string length and 
R is the curvature radius of the AdS space. In this limit of 
small string coupling and large curvature radius, classical 
gravity is a good approximation of the string theory. 

The background metric corresponding to the SYM the- 
ory at finite temperature is 



ds 2 



R 2 



du 2 



h(u) 



h(u)dt + u 2 (dx 



i\2 



G^dX^dX" 



with h(u) = u 2 — ujju 2 and where Uh = irT is a black- 
hole horizon. The corresponding Hawking temperature T 
is the temperature of the SYM theory. The coordinate 
in the fifth dimension u = r/R 2 has the dimension of 
momentum and the SYM theory lives on the boundary at 
u = oo. 

The heavy quark, in the fundamental representation, 
whose energy loss we want to compute, leaves on a brane 
at u = u m = 2irM/VX — ► oo with a string attached to 
it, hanging down to the horizon. Points on the string can 
be identified to quantum fluctuations in the heavy quark 
wave function with virtuality ~ u. Indeed, the quantum 
dynamics in the SYM theory is mapped onto classical dy- 
namics in the 5th dimension. More precisely, the string 
dynamics is given by the Nambu-Goto action 



5' 



2nR 2 



J drda^J- det g ab (17) 



where r and a are the worldsheet coordinates and 

g ab = G flv (d a X^(d b X l/ ) 
is the metric induced on the worldsheet. 



(18) 



3.2 The trailing string picture 

Let us assume that the quark moves along the x direc- 
tion and parametrize the space-time coordinates X^ = 
(t, x, y, z, u) — (r, x(t, a), 0, 0, a). Then 



S 



Vx 

2tt 



dtdu \ 1 — 



h(u) 



u 2 h(u)x' [ 



(19) 



with x = dtx(t,u) and x' = d u x(t,u). The classical equa- 
tion of motion d a 8C/ SdaX 1 — gives 



d ( u 2 h(u)x' 

du 



u d 
h(ujdi 







(20) 



From the solution of this equation, one gets the rate at 
which the energy flows down the string: 



dE 
~dt 



SC 

Wx 1 



VXR 2 , 



u h(u)xx' 



(21) 



which is identified to the heavy-quark energy loss. 

In Ref. 8,9 , the authors imagined using an external 
force to pull the quark at a constant velocity v. Writting 
x(t, u) = x + vt + F(u) they obtained (see Fig(2]) 



F(u) 



1 

2u h 



tan 



cot" 



and the corresponding rate of energy loss 



dE 
~dt 



VA 
2tt 



2 9 



(22) 



(23) 



(16) 3.3 Introducing the saturation scale 



Investigating this picture in more details, one finds that 
there is a special point on the string, at u = y^Uh (10] - 
In terms of quantum fluctuations in the heavy quark wave 
function, this is a special momentum scale that we shall 
denote Q s for reasons we explain now. 

First, the part of the string below y^Uh is not causally 
connected with the part of the string above: this point cor- 
responds to a black hole horizon in the rest frame of the 
string. Points on the string with u < y/jUh do not know 
about the heavy quark. Second, the energy density around 
the quark is unchanged up to distances l/(y/juh) : fluctu- 
ations with a virtuality higher than , localized within 
this circle, do not feel the plasma. In the fifth dimension, 
this means that points on the string with u > ^juh do 
not know about the plasma. 

Therefore, the part of the string above y/jUh corre- 
sponds to highly virtual fluctuations still part of the heavy 
quark wave function while the part of the string below 
y/jUh corresponds to longer-lived fluctuations which be- 
came emitted radiation. By analogy with the weak cou- 
pling picture, one is led to call this momentum scale the 
saturation scale 

Q s = Vju h . (24) 



4 



Cyrille Marquet: Heavy-quark energy loss in pQCD and SYM plasmas 




Fig. 2. The trailing string solution x(t,u) = xq +vt + F («). 
The part of the string above u = v / 7 1i '> is genuinely part of 
the heavy quark and the part of the string below ^fUh is 
emitted radiation. At u= ^yuh, the overlap in the longitudinal 
direction between the trailing string and a straight string is 



Note that with this understanding, it is straightforward to 
interpret the limiting velocity phenomenon M > V^7 T. 
The maximum value of 7 is reached when Q s crosses the 
heavy-quark brane u = u m . 

We now understand that quantum fluctuations in the 
heavy quark wave function become emitted radiation if 
k±_ < Q s - The maximum energy to a fluctuation can have 
can be estimated from the minimum longitudinal extend 
its need to overlap with the heavy quark. At u = Q s , which 
is where the highest u) will be reached, the overlap in the 
longitudinal direction between the trailing string and a 
straight string (which would be the solution in the vac- 
uum) is l/"/Q s which means u> < jQ s - These conditions 
on k± and uj are the same than in the pQCD case when 
written in terms of Q Sl therefore it is no surprise that the 
energy loss at strong coupling can be written (with v = 1 
at high energies) 



dE 
~dt 



(25) 



The only differences are that the probability for the fluc- 
tuation is now \/\ instead of a s 7V c , and of course Q s is 
different. This picture does not allow to determine prefac- 
tors, but it gives correct parametric dependences. 

Finally, one obtains the following relation between the 
coherence time of the dominant radiated fluctuations t c = 
j/Qs and the saturation scale Q s : 



Qs = t c T 2 



(26) 



This is in agreement with the scale found in [IT] (see 
also [12]) which separates weak from strong scattering in 
the deep inelastic structure functions calculations. Inter- 
estingly enough, \/Q s is also the screening length found 
in [13] which determines whether a quarkonium in the 
strongly-coupled SYM plasma is bounded (d < l/Q s ) or 
has dissociated (d > l/Q s ), with d the quark-antiquark 
separation. 



3.4 The case of finite-extend matter 

Addressing the case of finite-extend matter is important 
for phenomenology. An exact calculation is challenging, 
because if would require to extend the AdS/CFT corre- 
spondence to introduce the plasma length L in the metric. 
Instead, using the picture developed above, one can infer 
the result. The discussion is similar to the one in pQCD: if 
L > t c = y/j/T, the results are that of the infinite matter 
case and if L < t c , the energy of the fluctuations which 
dominate the energy loss is w = LQ 2 S and one has 



dE /T ^2 
~dt * Q s 



(27) 



with 



Qs = LT 2 . (28) 

Let us describe a calculation which further motivates these 
results [10] . A brief constant acceleration a to the desired 
speed v mimics the creation of a bare quark-antiquark 
pair. The accelerating string solution [14] reveals that u = 
a is a black hole horizon |15| . similarly to the Unruh ef- 
fect in general relativity. Highly virtual fluctuations cor- 
responding to u > a are part of the heavy quark while 
the ones corresponding to u < a are radiated due to the 
acceleration. This insures that if a > T, which we assume, 
the medium is not felt during the creation process. 

Then stopping the acceleration triggers the building 
of the wavefunction, essentially the separation at u = a 
decreases as j/t. The bare quark is turning into a dressed 
quark while interacting with the medium. The key issue 
is to understand the time it takes for the heavy quark to 
build the fluctuations which will be freed and those that 
dominate the energy loss. While in pQCD this was easily 
evaluated because of the picture of multiple scattering lo- 
cally giving transverse momentum to the radiated gluons, 
this is more subtle at strong coupling. The findings in [10 
support the results (|27M28|) . Once again this does not allow 
to determine prefactors, but the L and T dependences are 
robust. 



3.5 Results for -p^ broadening 

Let us give our related results for the p±— broadening of 
the heavy quark at strong coupling. Our picture gives 



dt dt 



(29) 



where t — t c in the infinite matter case and t — L in 
the finite matter case. Note that because vA ^> 1, ra- 
diative p±— broadening is dominant at strong coupling, 
as opposed to weak coupling. Also contrary to the pQCD 
case, the results are different for infinite and finite matter. 
In the former case, one obtains 



dpi 

dt 



A 7 T A 



(30) 



Contrary to the energy- loss formula (f2"3"]) . this result is 
non trivial to get with a direct calculation [16017] . but the 
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7 and T dependence come in a straightforward manner 
within our picture. In the finite matter case one gets 



at 



(31) 



Infinite and finite matter results are different at strong 
coupling because p? broadening is not a local phenonemon 
like in pQCD: there is no picture of local scatterings. At 
strong coupling, the transport coefficient q is not relevant, 
it does not control the energy loss or broadening. 



4 Conclusions 

Heavy-quark energy loss and p± — broadening have identi- 
cal parametric form when propagating through a weakly- 
coupled QCD plasma or a strongly-coupled SYM plasma, 
when written in terms of the saturation momentum Q s : 



dE 



dt 



a s N c 



dp 2 A 



dt 



a s N c 

Vx 



dQl 
dt 



QCD 
SYM ■ 

(32) 

This is not surprising as Q s , in both cases, is defined as the 
scale with controls which heavy quark fluctuations become 
emitted radiation. While for a pQCD plasma it is known 
that such a scale can be singled out, the fact that one could 
identify one also at strong coupling is the main message 
of this work. Then the different couplings a s N c and Va 
simply reflect the fluctuation probabilities in the heavy 
quark wave function. The saturation scale Q s is given by 

— for infinite-extend matter or L > t c = j/Q s 



Ql = (?7) 2/3 in PQCD and Q 2 , = T 2 7 in SYM 



— for finite-extend matter with L < t c 



(33) 



Ql = qL in pQCD and Q 2 S = T*L 2 in SYM . (34) 

We recall that q ~ a s T 3 in pQCD. 

On the right-hand side on the p±— broadening formula, 
t = t c in the infinite matter case and t = L in the fi- 
nite matter case. The pQCD result dp\jdt = a s N c q is 
for radiative p±— broadening and is actually subdominant 
(a s <C 1), p±— broadening is multiple scattering domi- 
nated in pQCD dp\jdt = q. In SYM, because A > 1, 
p±— broadening is radiation dominated. Hence this is a 
non local phenomenon and the transport coefficient q is 
not a relevant quantity at strong coupling. However if in 
the pQCD energy loss problem the soft scales in the pro- 
cess are subject to a strong effective coupling, q will be 
enhanced by that dynamics }18j . 

Finally, we obtained the first estimate of the plasma 
length dependence at strong coupling of the heavy-quark 
energy loss and broadening. The conclusion is that the 
L dependence in much stronger at strong coupling, (AE ~ 
L 3 ,Apl ~ L 2 ) in SYM while (AE ~ L 2 ,Ap\ ~ L) in 
pQCD. We also note the results discussed in [T5], that the 
energy loss at strong coupling will behave as AE ~ L, if 
the relevant dynamics is that of early times, before ther- 
malization, while in this work we considered a thermalizcd 
medium. 
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